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Optimal control theory implementations for an efficient population transfer and creation of a 
maximum coherence in three-level system are considered. We demonstrate that the half-STIRAP 
(stimulated Raman adiabatic passage) scheme for creation of the msiximum Raman coherence is the 
optimal solution according to the optimal control theory. We also present a comparative study of 
several implementations of optimal control theory applied to the complete population transfer and 
creation of the maximum coherence. Performance of the conjugate gradient method, the Zhu-Rabitz 
method and the Krotov method has been analyzed. 
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. ' Efficient and selective transfer of population in quantum systems provides an essential tool for a variety of applica- 
, tions in physical chemistry, laser spectroscopy, quantum optics, and quantum information processing [1- 5]. Some of 
the population transfer techniques make use of the Rabi oscillations vifhen the transfer efficiency is controlled by pulse 
area of the external fields [6] . Sensitivity of the final population distribution in the system to the field parameters is 
sometimes considered as a dravifback of these methods. Techniques utilizing adiabatic passage solution in the system 
dynamics are substantially more robust against moderate variations in the interaction parameters. Stimulated Raman 
^ ' adiabatic passage (STIRAP) is one example of such methods 7, 8]. 

d . In a three-level A system, STIRAP provides a robust scheme for transferring population by using two strong pulses, 
^ ' called the pump and the Stokes f?, ^ . When the pump-Stokes pulse sequence is arranged in a counterintuitive order 
(the Stokes pulse precedes the pump pulse) the system dynamics takes place within a so-called dark state (a coherent 
superposition of the initial and the target states). The adiabatic change of the field amplitudes guaranties nearly 
100% efficiency of population transfer to the target state while the intermediate state population is negligibly small 
during whole time evolution. There are a few modifications of the STIRAP scheme for more complex configurations 
of quantum systems (?- 15|. 
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\^ • Search for an efficient scheme of population transfer is also one of the major goals of optimal control theory (OCT), 

, a powerful and sometimes mathematically very sophisticated techniqu e | 1 6i-[20l. The OCT algorithm has been applied 

' successfully to a broad variety of physical and chemical systems 2l|-|30| . It intrinsically accounts for the quantum- 

, mechanical interference between many pathways connecting initial and target state in a quantum system. There were 

T-H ■ some interesting arguments in trying to find a link between OCT and STIRAP. Despite initial negative prognosis 13 1|. 

, it was demonstrated recently that STIRAP type solution can emerge automatically from global OCT method [32|. 

h>. ' Besides, STIRAP scheme has been obtained successfully using a local optimal control theory [1^ 



In this work we review a performance of the optimal control theory algorithms jH, 0, 0, l33l-36[ applied for the 



^ . population transfer and explore their potential to maximize the Raman coherence in the three-level A system. Our 
motivation for maximizing coherence is related to the recent developments in coherent anti-Stokes Raman (CARS) 
microscopy and remote detection using intense femto-second laser setups 37-4^. 



The paper is organized as follows. In section II, basic equations of the optimal control theory are developed using 
the method of variational calculus. Field equations are derived using the penalty on the energy of the control field. 
A time dependent penalty function is used which ensures experimentally feasible profile of the laser pulses. A second 
penalty function is introduced to minimize the population of the excited intermediate state throughout time evolution 
of the system. In section III, we apply the OCT formalism to a three-level A system and analyze solutions of the 
OCT equations for population transfer and a maximum Raman coherence applying different optimization strategies: 
the conjugate gradient method 2^ 3^1, the Zhu-Rabitz method [3] and the Krotov method 17 -2^. Section IV is 



the conclusion. 
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GENERAL EQUATIONS OF OPTIMAL CONTROL THEORY 



The OCT is based on the definition of a cost functional K which must have an optimal value when the desired 
transformation of the wave function is successfully achieved by the control laser field e{t). An optimal solution requires 
that the system wave function at a final time, 1-0(2^)), should be as close as possible to the target wave function, |0(T)). 
That is the overlap |(?A(T)|0(r))p is maximal at final time T. 

In order to derive control equation for the field and obtain realistic field amplitude we minimize the energy fluency of 
the field. Another requirement is that the population of the excited intermediate states has to be minimal throughout 
the transfer process 17, 3 32|, that can be done by defining a projection operator P = \4'intit)){4>intit)\ = J2k 
where is the eigenket of the unwanted intermediate states. There is an additional constraint that should be 
taken into account when we maximize or minimize the cost functional K: the wave function IV'Ci)) raust satisfy the 
Schrodinger equation. 

These requirements lead to a complete cost functional of the form 



K = |(0(r)|0(T))|' - a{t) / dt[eit) - 
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The factor a{t) — aQ/s{t) is a time dependent penalty function that determines the shape function s(t), aa is a 
constant which should be determined due to the significance of the field energy value. The main purpose of s(t) is 
to turn the field on and off smoothly to ensure feasibility of experimental implementation of the optimal laser pulse. 
er{t) denotes a reference field and /3 is a penalty parameter for the total population of the intermediate state (or 
state manifold). The function |x(^)) can be regarded as a Lagrange multiplier introduced to assure satisfaction of the 
Schrodinger equation. 

Each of the terms in Eq. ([1]) depends explicitly or implicitly on the unknown driving field. The goal is to determine 
an optimal field e{t) for which the cost functional has an extremum. Taking variations of the cost functional with 
respect to |x(^))> IV'(O) ^-^d e(0 we find the following set of equations 



ih-m)) = H\m), 



||xW) = -Ji?|xW) + WW> 



(2) 
(3) 



e(t) = er{t) + 



a(t)h 



lm(x{t) 



dH 



de{t) 



(4) 



Variation of the cost functional with respect to \^-){T)) gives the initial condition for the Lagrange multiplier 

\x{T)) = {^{Tmmm) ■ (5) 

Equations ([2]) and ([3]) determine the time evolution of the wave function and Lagrange multiplier which are used 
in Eq. dU to find the optimal field, e(t), maximizing the overlap |(V'(T)|0(r))p in Eq. ([T]). 



OPTIMAL CONTROL OF THREE-LEVEL A SYSTEM DYNAMICS 



General equations 



To demonstrate an application of the general optimization formalism outlined in the previous section and to test the 
performance of various implementations of the OCT we consider two optimal control problems. First, we reexamine 
the population transfer from the initially occupied level |1) to the level |3). In the second problem, we utilize the 
OCT to create a maximum Raman coherence, in other words, the 50/50 coherent superposition of states |1) and |3). 

First, we consider the population transfer in a generic three-level A system exited by a pump-Stokes pulse sequence, 
see Fig. [TJ We address the so-called nonimpulsive excitation when the pump pulse interacts with |1) — |2) transition 
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while the Stokes pulse controls the transition |2) — |3). We assume that all relaxation times in the system are much 
longer than the pulse duration so that the dynamics of an arbitrary wave function 



m)) ^ ai{t)\l) + a2{t)\2) + a3{t)\3) . 



(6) 



where ai{t) is the probability amplitude to be in state \i), is governed by the Schrodinger equation with the Hamiltonian 
of the form 



H 



El 
-M2iep(t) 




-Mi2ep(t) 

E2 





E3 



(7) 



here Ei is the energy of the i = 1, 2, 3 state, ^^12,23 £^re the dipole moments, ep,s(0 the pump and Stokes fields. 

Now we are ready to solve the optimization problem by applying a general algorithm outlined in the previous 
section. Varying the cost functional with respect to the Lagrange multiplier vector bi{t) , the probability amplitudes 
vector ai{t), and the external field €p^s{t), Eqs. ([U, ^ and ^ takes the following form: 



(8) 



dt 



'-H^jbj{t)+ I3a2{t)5a, 



(9) 



2a{t) 



Im 



bl{t)a2{t) + bl{t)ai{t) 



(10) 



At) 



2a{t) 



Im 



6*(t)a3(t) + fo3W«2(t) 



(11) 



where t^p{t) and e^(t) are the pump and Stokes reference fields. 

The set of Eqs. ([8|) - ([TT|) provides optimal control fields which are based on the solution of the Schrodingrer 
equation with the Hamiltonian in Eq. ([T]). In many cases, we consider the interaction of the three- level system with 
the external field as a resonant excitation and invoke the Rotating Wave Approximation (RWA) which partially 
simplifies the solution of the problem and sometimes facilitates understanding of underlaying physical mechanisms. 

In the interaction representation the Hamiltonian can be written in the form 
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where W^p,s(i) = ilp,s(i)(l + e~^*'^^'^*)> ^p,s(0 = lJ'i2,23E^s{t)/h are the pump and Stokes Rabi frequencies, E^g{t) 
are the envelope of the pump and Stokes pulses, Ap^s — ^p,s ~ 1^21,23 are the single photon detunings of the pump 
and Stokes central frequency ujp^s from the respective transition frequency a;2i.23- 

Using the RWA we can replace Wp g{t) in Eq. p2|) by the corresponding Rabi frequencies, that is we neglect the 
rapidly oscillating terms. As a result we rewrite the equations for the optimal fields, Eq. ([T0|) . in terms of the 

Rabi frequency envelopes of the pump and Stokes fields 



1 



2a(t) 



Im 



b{{t)a2{t) + b*{t)ai{t) 



(13) 
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bl{t)a^{t) + bl{t)a2{t) 



(14) 



where fip(t), i^5(i) are the reference Rabi frequency of the pump and Stokes pulses. 
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FIG. 1: Schematic of a three-level system. 



Complete population transfer 



To examine OCT implementation we consider three different optimization methods 
method 



33l . |43|, the Zhu-Rabitz method [16| and the Krotov method [17|, |l 



the conjugate gradient 
A detailed description of the 
numerical schemes is given in the Appendix. We choose the Gaussian form for the pump and Stokes pulse envelopes 
as an initial guess 



Qp g[t) = rig exp 



ZTq 



(15) 



The reference envelope 



and our target time is equal to T — 10 in the units normalized by the time duration tq. 
rip g{t) = is used unless otherwise stated. 

The goal of the first problem at hands is twofold: first, to design the shape and the sequence of the pump and 
Stokes pulses providing a complete population transfer to state |3); second to suppress population of the intermediate 
state during excitation process by applying a penalty on the state population. By doing this we also would like to 
determine the efficiency of the methods mentioned above. For simplicity we restrict our consideration to the exact 
resonance conditions, Ap — As — 0. 



The conjugate gradient method. The results obtained using the conjugate gradient method [29|, |33|, |43[ are 
shown in Fig. [2l The set of plots (a)-(c) shows results obtained without the penalty on the intermediate state 
population, (3 = 0. Plots (d)-(f) show results obtained when the penalty is imposed on the population of state 
|2), /3 ^ 0. The population of the states |1), |2), and |3) as a function of time is shown in plots (a) and (d). 
Complete population transfer from initial level |1) to the target level |3) is achieved at the final time T independently 
of the penalty parameter. The amount of population in the intermediate state is considerably reduced for a field 
obtained with the state-dependent constraint (Fig.[2Ud)) in contrast to that resulting from unconstrained optimization 
(Fig. (Ha)). The optimized Rabi frequencies (obtained after 1000 iterations) are shown in Fig. [2] (b) and (e). With 
a proper choice of the penalty parameters and /?, we remove almost all the intuitive solutions and obtain the 
STIRAP solution, a counterintuitive pulse sequence when the Stokes pulse precedes the pump. Plots (c) and (f) of 
Fig. [5] show the convergence behavior of the optimized transition probability defined as V = |(V'(r)|(/)(T))p, and the 
cost functional, K , versus the number of iteration steps. It is seen that the transition probability reaches nearly 100%. 



The Zhu-Rabitz method. Figure [3] shows the dynamics of population transfer in the three-level A system optimized 
by using the Zhu-Rabitz method [l6|. The same order of illustrations and legends is kept in Fig. |3]as in the case of 
Fig.H 

There is a complete transfer of the population from the initial state to the target state at final time T, Fig. [3] (a) 
and (d). However, we observe two different mechanisms of population transfer. Optimization without a penalty on 
the intermediate state population (left panel of Fig. [3]) produces an intuitive pulse sequence, first the pump, then the 
Stokes pulse is applied to the system (Fig. [3Jb)). That pulse arrangement results in a sequential population transfer 
from the state |1) to the state |2) and then to the target state |3) as it is demonstrated in Fig. [3]^ a), with almost 80% 
of population resided in the state |2) at the central time. 
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FIG. 2: Population transfer in a A system using conjugate gradient metliod with penalty on population of the state |2) (right 
panel) and without penalty, /3 = (left panel), (a), (d) The population of state |1) - red line, |2) - blue line and |3) - black line, 
(b), (e) The sequence of optimal pulses: Rabi frequency of the pump pulse - red line and Stokes pulse - blue line, (c), (f) The 
convergence behavior of the optimized transition probability - red line, penalty on the field energy - blue line, penalty on the 
intermediate state population - pink line and final optimized cost functional - black line versus the number of iteration steps. 



Applying a penalty on the second state population changes the mechanism dramatically (see right panel of Fig. [3]). 
Now we obtained almost hundred times more intense pulses, compare Fig. [3] (b) and (e). More importantly, we 
automatically obtain a counterintuitive pulse sequence when the Stokes pulse precedes the pump pulse. Therefore, 
it is clear that the OCT algorithm finds a well known STIRAP solution which has a highly pronounced signature of 
suppressed intermediate state population, see Fig. |3] (d). In that case population transfer takes place through the 
dark state consisting ideally only of the initial and the target state probability amplitudes and has no projection to 
the state |2) 0. 

The value of the optimized cost functional in Fig. [3] (f) is a bit less than that of the cost functional value obtained 
by using the conjugate gradient method. Fig. [5] (f). This is due to a slightly larger energy fluency of the fields as it is 
clear from the field amplitudes in the plot (e) in Fig. |3] comparing to that in Fig. [5] 

The Krotov method. The optimized results of population transfer using the Krotov method 17|, |l8| are shown 
in Figs, m and [S] The same set of parameters is used in the optimization process as in the previous two cases. Figure[S] 
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illustrates the results obtained using the Krotov method when f2pg(f) ^ Population of the states, optimized 
Rabi frequencies, transition probability and cost functional are shown in plots (a)-(c) and (d)-(f) of Figs. |4] and [5j 

As in the previous cases, we produce complete population transfer to the target state. However, using Krotov 
method 17], |l8| with fip g{t) — we were not able to find a set of penalty parameters which provides considerable 
suppression of the second state population, see Fig. |3](d). At most we are able to reduce the second state population 
to only about 40%. Optimized Rabi frequencies are arranged in the intuitive manner even for the reasonably large 
areas of the pulses, Fig. |4] (b) and (e). 

The situation changes considerably when we use the reference fields, il^p g{t) ^ 0, Fig. [S] In the manner similar to 
previous two cases the STIRAP solution emerges from the optimization procedure when we impose a penalty on the 
intermediate state population, see right panel in Fig[5l the Stokes pulse precedes the pump Fig[5](e) and a detrimental 
population of the state |2) is suppressed almost to zero, Fig[S](d). 

Table U shows a comparison of transition probability 7^, optimized cost functional K and the maximum population, 
Qt2 = 1^2 which resides in the intermediate state at central times for all the methods. Note that the value of the cost 
functional obtained using the Krotov method for fl^p g{t) 7^ is larger than corresponding value in other methods. 
However, the computation cost of the conjugate gradient method is larger than that of the Zhu-Rabitz method and 
the Krotov method. 
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FIG. 3: Population transfer in a A system using Zhu-Rabitz method; description of plots same as in Fig. [2] 
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Optimal control to maximize coherence 



In this section we apply the OCT to create a maximum coherence \g3i\ — \a^o,i\ between the levels |1) and |3) at 
final time T. We use the same parameters and the initial guess function for the fields envelopes as in the complete 
population transfer section. Figure |6] illustrates the results obtained using the conjugate gradient method (a)-(c), the 
Zhu-Rabitz method (d)-(f) and the Krotov method with ftp g{t) ^ (g)-(i), respectively. This time we restrict our 
consideration to the case when a penalty on the intermediate state population is applied. 

Figure [6] (a), (d), (g) shows the time evolution of the population in the three- level system excited by the optimized 
pulse sequence. Dynamics of the population presented in the Fig. [6] is almost identical for the all three versions of 
the OCT: at the target time we obtain a maximum coherence |£>3i(T)| 1/2, that is we create a 50/50 coherent 
superposition of states |1) and |3). Population of the intermediate state |2) is almost negligible during the excitation 
process. 

The Rabi frequency of the pump and Stokes pulses obtained from the implementations of the OCT is shown in 
Fig[6](b), (e) and (h), respectively using the conjugate gradient method 2^ H, 43|, the Zhu-Rabitz method [3| and 
the Krotov method with fip ^ (t) 7^ [19| . All three methods give a similar solution for the optimal pulse sequence. 




FIG. 4: Population transfer in a A system using Krotov method: description of plots same as in Fig. (2] 
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The obtained pulse sequence and corresponding population dynamics allow us to conclude that the OCT finds the 
so-called half-STIRAP (also sometimes referred as fractional STIRAP) scheme as the optimal pulse sequence to create 
maximum coherence in three-level system. The same way as in the STIRAP scheme, the Stokes pulse is turned on 
first but both pulses are turned off simultaneously at the later time, see Fig[n](b), (e) and (h). The mechanism of the 
solution can be explained using the dressed state basis as follows. It is easy to find the eigenvalues and eigenvectors 
(dressed states) of the system Hamiltonian, Eq. in the RWA. The important dressed state with energy Aq = 

has the form 

\co{t)) - [ns{t)/^nut) + ni{t), o, -npit)/^nj,it) + niit)}. (i6) 

That dressed state has zero projection on the pure state |2) while probability amplitudes to be in the state |1) and 
|3) are controllable by the ratio of the pump and Stokes Rabi frequencies. Analyzing Eq. (IT6|) we can reproduce the 
population dynamics for the half-STIRAP scheme: at time t = the pump Rabi frequency ^lp{t = 0) = and the 
Stokes Rabi frequency ^si^) ^ therefore the dressed state \co{t = 0)) correlates with the initial state during the 
turn off stage of the pulse excitation ratio between pump and Stokes Rabi frequencies is equal to one, that results in 
creation of state \co{t = T)) = 1/a/2, 0, — 1/\/2, the state of maximum coherence between states |1) and |3). Whole 
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FIG. 5: Population transfer in a A system using the Krotov method for Sip g{t) 0; description of plots same as in Fig. [21 
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TABLE I: Comparison of results obtained using difTerent implementations of the OCT for the population transfer dynamics in 
the three-level A system. 



Method 


Without state-dependent penalty 


With state-dependent penalty 


ao 


/3 


V 


K 


£122 


Qo 


/3 


P 


K 


£122 


Conjugate gradient 


0.01 





0.999 


0.836 


0.48 


0.00005 


1.0 


0.998 


0.930 


0.04 


Zhu-Rabitz 


0.01 





0.999 


0.882 


0.74 


0.0005 


1.8 


0.998 


0.805 


0.053 


Krotov {Q.pA^') = 0) 


0.01 





0.999 


0.883 


0.74 


0.005 


0.2 


0.998 


0.805 


0.47 


Krotov (np,s(i) / 0) 


1.0 





0.999 


0.999 


0.46 


0.05 


0.2 


0.999 


0.996 


0.024 



dynamics of the system takes place in one dressed states and that is possible only in the limit of the adiabatic regime 
when nonadibatic coupling to other two dressed states is negligible. According to our observation the OCT finds the 
adiabatic mechanism as an optimal solution of the control problem. 

To compare the performance of the different implementations the optimized transition probability 7^, the final 
optimized cost functional K, the penalty on the field energy and the penalty on the second state population as a 
function of number of iterations are shown in Fig|6] (c), (f) and (i). Best results for the maximal coherence obtained 
using the discussed methods are also summarized in Table HT] 




FIG. 6: Maximizing coherence in a three-level system using conjugate gradient method (a)-(c), Zhu-Rabitz method (d)-(f) 
and Krotov method (g)-(i). (a), (d), (g) The population of state |1) - red line, state |2) - blue line, state |3) - black line and 
coherence £131 - pink line, (b), (e), (h) The sequence of optimal pulses: Rabi frequency of the pump pulse - red line and the 
Stokes pulse - blue line, (c), (f), (i) The convergence behavior of the optimized transition probability - red line, penalty on 
the field energy - blue line, penalty on the second state population - pink line and final optimized cost functional - black line 
versus the number of iterations. 
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TABLE II: Comparison of the results obtained using different implementation methods of the OCT to maximize coherence, 
|p3i(r)|, between levels |1) and |3) in the three-level A system; Qii{T) — \ai{T)\'^ is the population of state \i) , i = 1,2,3. 



Method 


QO 




£-11 


£122 


£133 


l£'3l| 


V 


K 


Conjugate gradient 


0.00025 


0.2 


0.521 


0.058 


0.478 


0.499 


0.999 


0.973 


Zhu-Rabitz 


0.0005 


1.8 


0.538 


0.030 


0.461 


0.498 


0.998 


0.905 


Krotov / 0) 


0.1 


0.2 


0.503 


0.011 


0.496 


0.499 


0.999 


0.997 



CONCLUSION 



In this paper, we have analyzed a performance of several implementation methods of the OCT to design optimal 
pulse sequences for a complete population transfer and the creation of a maximum coherence in a three-level A 
system. We have applied the conjugate gradient method 2^ 33, i^l, the Zhu-Rabitz method fli'l, and the Krotov 



method [17H20j to obtain the optimal solution. 

Using the conjugate gradient method it was demonstrated earlier [l^l that STIRAP type solution can emerge au- 
tomatically from the global OCT method. Now we have shown that the Zhu-Rabitz method [ig'] and the Krotov 



method [17h20| provide an additional conformation that counterintuitive pulse sequence is the optimal solution ac- 
cording to the OCT. It was demonstrated that the penalty on the population of the intermediate state in a three-level 
system is a crucial factor to obtain the optimal STIRAP scheme. 

We have also demonstrated that the half-STIRAP scheme is the optimal solution according to the optimal control 
theory for creation of a maximum coherence. We have shown that a creation of a maximally coherent superposition 
is one more optimization problem where the OCT finds the adiabatic mechanism as an optimal solution. 

Note that the related results using OCT to create predetermined coherent superposition has been published re- 
cently (29j . However, the target wave function used in [23w_as different in comparison to present studies: population 
of the initially populated state at the later time was zero 29]. As the result, the generalized STIRAP pulse sequence 
was found as the optimal solution using the conjugate gradient method (29| . In our case, the optimization procedure 
reveals the half-STIRAP pulse sequence as the optimal solution since we require that the half of the population should 
remain in the initially populated state. Comparing the results obtained using the above considered methods, it is 
clear that the methods can be applied to various problems on equal footing, since the value of optimized transition 
probability V is more than 99% for all the methods. However, the value of optimized cost functional K calculated 
using the conjugate gradient method is larger than corresponding values calculated using the Zhu-Rabitz iterative 
method and the Krotov method for ^"^p g{t) = and /3 7^ 0. Note that the computation cost of the conjugate gradient 
method is generally larger than that of the other two methods. 

One might observe that the pulses in optimal sequence obtained by optimizing the population transfer (Figs. 2-5) 
and the maximum coherence (Fig.6) have a relatively similar structure in all considered implementation schemes. 
However, the Zhu-Rabitz method [l6| provides the shorter and more intense pulses in comparison to other methods. 
There is also a bit more pronounced asymmetry of the pump-pulse shape obtained by the Krotov method. We believe 
that these differences are due to the variations in the numerical implementation procedure and the values of the 
penalty parameters. The conjugate-gradient method 2^ 33, 4^ incorporates one iteration to next iteration feedback 
from the control field while Krotov method 17-2^ incorporates one time step to next time step feedback of the 
control field, and Zhu-Rabitz method 1^ incorporates feedback from the control field in an entangled fashion [45| . 
mixing up both feedback techniques. The detailed description of the implementation procedure of these methods is 
given in the Appendices. 
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Conjugate gradient method 



For given target time T and the number of time steps N {ti — i x At, where i = 0,1, 2, N and NAt — T), the 
conjugate-gradient method 2^, 33, 43| involves the following steps to obtain the optimal field: 
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Step 1: Choose an initial electric field eo(<). 
Step 2: Set fc = 1 and e''{t) = eo{t). 

Step 3: Propagate \ip''{t = 0)) forward in time using e^-^^) according to Eq. 1^ to obtain \ip''{T)). 

Step 4: Evaluate the cost functional K'^ according to Eq. ([Ij. The last term in this equation is zero as \x{t)) is a 
solution of the time-dependent Schrodinger equation. 

Step 5: If fc > 2, compute AK'^ = — K^""^ and compare it with the convergence threshold, 7. If AK^ < 7, 
then stop the iteration and declare that the optimal pulse has been obtained. 

Step 6: Set {x''{T)\ = {<p{T)\{(l){T)\'ip'' (T)) . Propagate {x''{T)\ and \ip''{T)) backward in time using field e^{t) to 
obtain (x'^WI and {^''{0)). 

Step 7: The gradient g''{t) of the cost functional K'' defined in Eq. ([T|) with respect to the variation of e'^(t) at 
time t is given by 



de^{t) 



a{t)e'^{t)-lm{x\t) 



dH 



At. 



(17) 



The Polak-Ribiere-Polyak search direction is calculated using the Eq. ([T7| as 



where 

k — 2,3,..., d^{t) = g^{t). The function ('^ is the conjugate gradient update parameter and g'^{t)'^ is the transpose of 
g^{t). A line search is then performed along this direction to determine the maximum value of the cost functional. 
Step 8: The electric field for the next (fc + 1) iteration is taken as 

e^+i(i) = e^(i) + A-rf^(t), 

where A is determined by the line search, which makes e^'^^{t) to generate the maximum value of K. 
Step 9: Go back to Step 3 and repeat 3-8 until the required convergence has been achieved. 



Zhu-Rabitz method 

Zhu-Rabitz method 16] involves the following steps to find the optimal value of the control field: 
Step 1: Choose an initial electric field ei{t). Set fc = 1 and e^{t) = eo(t). 

Step 2: Propagate |'0'^(O)) forward in time using e'^{t) according to Eq. ([2|) to obtain \'4>^{T)). 

Step 3: Evaluate the cost functional K'^ according to Eq. If fc > 2, compute AK^ — — K^^^ and compare 
it with the convergence threshold, 7. If AK^ < 7, then stop the iteration and declare that the optimal pulse has been 
obtained. 

Step 4: Propagate Lagrange multiplier (x'=(r)| = (0(r)|(</)(T)|V'''(T)) backward in time using the new electric 
field e'=+i(t) defined by 



it) 



1 



ait) 



iMx'it) 



dH 



de^'+^t) 



according to Eq. ^ to obtain (x'''(0)|- 

Step 5: Propagate \il}^^^{{))) forward in time using the new field e'^''"^(t) given by 



(i) 



a{t) 



luv{x\t) 



dH 



de^+^{t) 



according to Eq. Q to obtain 1^;''+'^ (T)) . 

Step 6: Go back to Step 3 and repeat 3-5 until the required convergence has been achieved. 
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Krotov method 

The Krotov method [13, [3| employs a shghtly different iteration procedure: 

Step 1: Choose an initial electric field ei{t). 
Step 2: Set fc = 1 and e''{t) = e^{t). 

Step 3: Propagate \ip''{t = 0)) forward in time with the field e''{t) according to Eq. ^ to obtain \ip'^{T)). 

Step 4: Evaluate the cost functional K'^ according to Eq. If fc > 2, compute AK'^ = — K^~'^ and compare 
it with the convergence threshold, 7. If AK'^ < 7, then stop the iteration and declare that the optimal pulse has been 
obtained. 

Step 5: Set (x''(r)| = ((/>(T)|(0(r)|?/''' (T)>. Propagate (x''(r)| backward in time according to Eq. using field 
e'=(i) to obtain (x''(0)|. 

Step 6a: Set = 0)) and propagate forward in time according to the Schrodinger equation ([2|) with simul- 

taneous evaluation of the electric field e''~^^{t) at each time step. To calculate {x{t)\ use the old electric field, i.e. the 
electric field generated in the previous iteration step and calculate \ip'''^^{t)) with the new electric field. 

Step 6b: Compute the new electric field e'^'^^(t) from (x'^(t)| and \tp'^^^{t)) according to 

Step 7: Go back to Step 4 and repeat 4-6 until the required convergence has been achieved. 
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